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LIST OF SYMBOLS 



A 


2 

Coefficient of the x term in the equation of a 

line on a constant life fatigue diagram where 
minimum-stress is x and maximum stress is y. 

(R = Ax z + Bx + C - y) 


AA 


An assigned value of +1. or -1. 


AAA 


A stress used in the calculation of plastic strain. 


ABDIF 


The absolute value of DIF. 


ABM 


The absolute value of ASMAX or of ASMIN, as assigned 


ABMAX 


The absolute value of ASMAX. 


ABMEAN 


The absolute value of ASMEAN 


ABM IN 


The absolute value of ASMIN. 


ABR4 


The absolute value of R(4) . 


ABR7 


The absolute value of R(7) . 


ABS 


The name of a routine calling for the absolute value 
of a quantity. 


AKT 


Stress concentration factor, K^ . 


ASMAX 


The product (AKT) (STMAX) . 


ASMEAN 


The quantity (ASMAX + ASMIN) /2. 


ASMIN 


The product (AKT) (STMIN) . 


AVSGMN 


Average value of S-IGMIN over an interval. 


AVSGMX 


Average value of SIGMAX over an interval. 


B 


Coefficient of the x term. (see A) 


BBB 


A stress used in the calculation of plastic strain. 


C 


The constant. (see A) 


COFMAN 


Inverse of the Coff in-Manson slope. 



CYCINT 



The number of cycles in an interval. 



CYCLES 


The calculated number of cycles expected to be 
indicated on a constant life fatigue diagram for the 
applied combination of maximum and minimum stress. 


Cl 


The Residual Stress Relaxation Constant. (see ENEP) 


DAM 


Damage . 


DECK 


Decimal or real value of integer K after conversion. 


DEL 2 


A portion of a least - squares -method solution. 


DIF 


The difference between residual stress and equilibrium 
residual stress. (RES(I) - EQRES) 


D02 


A portion of a least-squares-method solution. 


DUMMY 


A variable used in the calculation of the number of 
cycles to be considered as an interval for relaxa- 
tion determination. 


ELMOD 


The elastic modulus. 


EN 


The number of cycles from the beginning of the relax- 
ation process to the end of the current interval. 


ENEP 


The number of cycles required for overload residual 
stress effect to return to within one-tenth of its 
original difference from equilibrium conditions. 

(N = C1/(ABM) E1 (ABMEAN) E2 ) 

ep 


ENN 


The number of applied cycles at a load level. 


ENNCYC 


The ratio of the number of applied cycles to the 
number of cycles to failure. (ENN/CYCLES) 


EPSD 


LCF strain intercept. 


EQRES 


Equilibrium residual stress. 


EX 


An exponential function depicting the relaxation 
of residual stress. 


EXP 


The name of a routine calling for the exponential 
value of a quantity. 


EXPO 


An exponent. The power of 10 which indicates the 
number of cycles to failure. 


El 


Residual stress relaxation exponent. 


E2 


Residual stress relaxation exponent. 
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FLOAT 


The name of a routine calling for integer- to-real 
conversion . 


I 


A variable subscript. 


IBLOCK 


The identifying number of a block, the blocks being 
numbered consecutively from 1 to NBLOCK. 


IFIX 


The name of a routine calling for real- to - integer 
conversion . 


IN 


The number of steps input to the range-pair counting 
subroutine . 


IPRINT 


Value controlling the WRITE statements. 


IRAIN 


A counter. 


IRPCM 


Value controlling entry into the range-pair counting 
subrout ine . 


I STEP 


The identifying step number, the steps being numbered 
from 1 to NLEVEL. 


ITYPE 


The identifying type number, the types being numbered 
from 1 to NTYPE. 


J 


A variable subscript. 


JA 


Value of +1 or 0, as assigned for branch determination 


JB 


Value of -1 or 0, as assigned for branch determination 


JJ 


An index variable. 


JJJ 


An index variable. 


JKL 


An index variable. 


K 


An index variable. 


KK 


An index variable. 


KPMAX 


The number of steps output from the range-pair 
counting subroutine. 


L 


An index variable. 


LMN 


An index variable. 


M 


An index variable. 


N 


An index variable with values of N=4-7 indicating 
the power of 10, and thus identifying a particular 
life cycle curve. 
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NBLOCK 


The total number of times to execute a block of 
loads . 


NDECK 


The number of data decks to be run sequentially. 


NFLAG 


An integer used as a counter. 


NFLAG2 


An integer used as a counter. 


NLEVEL 


The total number of steps, or levels, of loads in 
a block. 


NN 


A subscripted variable used to indicate which types 
of loads are experienced in which blocks. 


NTYPE 


The total number of different types. 


PLSTRA 


Plastic strain. 


R 


Residue term in damage calculation. 


RES 


Residual stress. 


RNCYC 


The number of cycles for a level after exitting 
the range-pair counting subroutine. 


SIGMAX 


Maximum stress. 


SIGMIN 


Min imum stress . 


STMAX 


Maximum applied stress. 


STMIN 


Minimum applied stress. 


SUMDEL 


Summation of damage for a flight. 


SUMENN 


Accumulated total of applied cycles. (Summation of 
ENN) 


SUMNC 


Accumulated cycle ratio. (Summation of ENN/CYCLES) 


SUMR 


Summation of R(N) , N=4,7. 


SUMRN 


Summation of NR(N) , N=4,7. 


SUMR 2 


Summation of R(N)", N=4,7. 


SUMR2N 


Summation of NR(N)^, N=4,7. 


SUMR 3 


Summation of R(N)^, N=4,7. 


SUMR4 


Summation of R(N)^, N=4,7. 


TITLE1 


Identification of the source of the S-N data. 
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TLL 



Tensile limit load. 



TM1 , TM2 

TTYS 

TYS 

T1,T2 

X 

Y 



Material type. 

One-fifth of tensile yield stress. 
Tensile yield stress. 

Test identifying information. 
Variable equivalent to SIGMIN. 
Variable equivalent to SIGMAX. 
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I . INTRODUCTION 



In an aircraft fatigue life monitoring program, deter- 
mination of the actual loading environment which produces 
the fatigue is a major problem. With both military trainer 
and fighter aircraft, which may be used for a variety of 
duties, regular or continuous load recording programs have 
to be considered mandatory. Such load recording programs 
provide the data to calculate the consumed fatigue life of 
operating aircraft and are used in selecting design spectra 
for future aircraft designs. 

To facilitate the calculations of fatigue life from 
such data, it is important to know not only the loading 
magnitudes but also the time sequence of the loading. The 
fatigue life of a local critical point varies due to 
residual stresses remaining after the application of a load 
causing local plasticity in tension or compression. For 
instance, a peak tensile load into the plastic flow region 
leaves a residual compressive stress, which lowers the 
local magnitude of the follo\ving tensile stresses, thus 
increasing the fatigue life. In a similar manner a peak 
compressive load into the plastic flow region, if it were 
to occur, would leave a residual tensile stress, which 
decreases the fatigue life. 

Currently, fatigue monitoring of naval aircraft is 
based on the total number of g readings recorded at four 
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selected levels by an exceedence level counting accelerometer 
Using micropressors in the near future it will be possible 
to record each maximum and minimum load level experienced 
by an aircraft in sequence. The data collected can be used 
to monitor the fatigue life of a structure via the deter- 
mination of damage accumulated at a point found to be 
critical in a structural test of a prototype. It should 
be noted that this calculation uses a theoretical model. 
Damage is not observable or measurable. 

The objective of this thesis is to use a computer program 
employing the range-pair counting technique on time sequence 
recorded maximum and minimum aircraft loads and a relaxation 
model to consider residual stresses, to calculate plastic, 
elastic, and total aircraft structural damage. 
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II. GENERAL DESCRIPTION OF THE COMPUTER PROGRAM 



The computer program used in this fatigue analysis is 
divided into four modules. Each module is clearly labeled 
in the program listing (page 93) . In module I the input 
data are read and assimilated in preparation for future 
calculations. The input parameter requirements are presented 
in detail beginning on page 18 , and six sample sets of 
input data are illustrated starting on page 40 . 

When the load sequence randomization technique is used 
to input loading data, it is a part of module I. The load 
sequence randomization technique uses the computer library 
subroutine RANDU to place 10 percent of the MIL spectrum 
A positive loads (Table I) in a random order. Each of the 
positive loads is paired with a minimum load of 11 percent 
limit load, or 1-g. During the randomization each load has 
an equal probability of selection. A counter restricts the 
number of times a value is selected to the number of 
occurrences of the particular load level in MIL spectrum A. 

Two variations on the randomization of the MIL spectrum 
A load levels are available. The first option permits 
randomization of the negative MIL spectrum A load levels 
using the same technique as in the case of the positive 
MIL spectrum A load levels. Each negative load level is 
randomly paired with a positive load level. Since there 
are not as many negative load levels in MIL spectrum A as 
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TABLE I 



FREQUENCY OF MANEUVER LOADS 

NUMBER OF TIMES PER THOUSAND HOURS THAT 
LOAD FACTOR IS EXPERIENCED 



PERCENT OF MAXIMUM (POSITIVE) FLIGHT MANEUVER LOAD 

SYMMETRICAL LIMIT LOAD FACTOR SPECTRUM A 



35 


17000 


45 


9500 


55 


6500 


65 


4500 


75 


2500 


85 


1360 


95 


440 


105 


150 


115 


40 


125 


16 

TOTAL 42006 



PERCENT OF MAXIMUM (NEGATIVE) FLIGHT MANEUVER LOAD 

SYMMETRICAL LIMIT LOAD FACTOR SPECTRUM A 



0 




500 


10 




200 


20 




100 


30 




60 


40 




35 


50 




30 


60 




25 


70 




20 


80 




15 


90 




10 


100 




5 


110 




3 



TOTAL TWI 
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there are positive load levels, the excess load levels are 
paired with 1-g loads. Mixing of the negative MIL spectrum 
A load levels and the 1-g load levels is accomplished using 
the computer library subroutine RANDU. For this ordering 
there is an 80 percent chance that each negative load selected 
will be a 1-g load. A counter ensures the proper numbers of 
each load type are included. Eighty percent probability of 
selection of a 1-g load was used to spread the smaller 
number of negative MIL spectrum A loads throughout the 
sequence, and also to speed computer operation by not having 
the matrix of negative MIL spectrum A loads addressed so 
often, after the counter indicated the matrix elements had 
all been used. 

The second option for randomizing loads provides for the 
inclusion of a ground cycle between each flight. The ground 
cycle loading is taken as -2500 psi, or -8 percent limit 
load. Eight percent is the level of the ground cycle load 
used by the A7 manufacturer. There are 4201 positive 
loading events in 100 hours of MIL spectrum A flight time. 

Each flight is arbitrarily taken as one* hour in length, so 
a ground load is paired .with every forty-second positive 
loading event. 

In module II the local stresses and strains are deter- 
mined from the input data. If the local peak stress is in 
the plastic region, the specimen is assumed to unload 
elastically, leaving a local residual compressive stress 
in the material. In the analysis developed by Potter [Ref. 6], 
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the transient portion of the residual stress so produced 
relaxes toward zero, or an equilibrium residual stress. 

After the stresses and strains are calculated, the local 
stress cycles are counted in module III using the range- 
pair counting technique. In module IV the fatigue damage 
is calculated. Damage is determined separately for elastic 
and plastic strain events. The elastic and plastic damage 
for each flight is listed as output, as well as the total 
cumulative damage from any previous flights. 

A. INPUT DATA REQUIREMENTS (MODULE I) 

Samples of data are listed starting on page Lo; format- 
ting is described beginning on page 93 . Both pages should 
be consulted while reading the following descriptions. 

1 . Data Card 1 

Only one of these cards is required or permitted 
for each program run. 

NDECK = the number of data decks to be run sequen- 
tially. The input on this card will be common to all data 
decks run. It is not necessary for the different data 
decks to have any parameters in common. 

IPRINT = the value controlling the write statements. 
Output available during the process of the analysis 
includes : 

a. Maximum and minimum applied stress of each cycle 
and the local stress response throughout the spectrum. Also 
printed out is the residual stress, equilibrium stress, 
applied cycles, and the equilibrium period for each cycle. 
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b. The elastic local stress history as input into 
the range-pair subroutine and the resulting range-paired 
spectrum . 

c. The maximum plastic strain occurrence during the 
spectrum and the damage associated with each strain reversal. 

d. The accumulated damage associated with the plastic 

strains . 

e. The range-paired elastic stress spectrum and the 
damage associated with each level for block loading, or for 
each cycle for RANDU generated loadings. 

f. The accumulated damage associated with the cur- 
rent block of loading, including the plastic strain damage, 
and the total damage since the initiation of cycling. 

The value assigned to IPRINT controls which pieces 
of data are included in the printout. 

If IPRINT = 1, all six items listed above are 
printed for each flight or block of loads. 

If IPRINT = 2, all items listed above except b, 
are printed. 

If IPRINT = 3, only items d and f listed above are 

printed. 

An example of each output option is given beginning 
on page . When the number of loads is more than fifty, 
the output is quite voluminous, and a paper usage default 
may interrupt program operation. With large numbers of 
loads, IPRINT should be set equal to 3. 



IRPCM = the value controlling entry into the range- 



pair counting subroutine. 



If 


IRPCM = 1, 


the 


range-pair 


counting 


method is used. 


If 


IRPCM = 2, 


the 


range-pair 


counting 


method is 



skipped . 

Each data deck run sequentially must contain the 
following cards. The number of data decks was indicated by 
the value assigned NDECK. 

2 . Card 1 

Information on this card is descriptive, alphanumeric 
data used only as a heading test identification. 

3 . Card 2 

The name of the material type and four material con- 
stants are listed on this card. The material type is an 
alphanumeric entry. The tensile yield stress, LCF strain 
intercept, inverse of the Coff in-Manson slope, and the 
modulus of elasticity are real number entries. The LCF 
strain intercept is the numerical value of the ordinate 
intercept in the log-log Coff in-Manson plot of plastic 
strain range vs. cycles to failure. Units of tensile yield 
stress and modulus of elasticity are ksi. 

4 . Cards 3 Through 6 

Each card contains the three coefficients of a second 
order leas t- squares curve fit for S-N data obtained from a 
Goodman diagram published in MIL-HDBK-5A. The curves are 
fit for lives of 10^, 10^, 10^, and 10 ^ cycles. Space is 
provided on each card for an alphanumeric indication of 
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the data source. Only the data source printed on card 6 
will be retained in computer storage for subsequent printout 
Columns 72 through 80 are used for information only and are 
not read by the computer. 

5 . Card 7 

Constants to be used in the equilibrium period cal- 
culations are listed on this card. 

6 . Card 8 

The stress concentration factor is the only entry 
on this card. 

7 . Card 9 

The first entry on this card is the number of blocks 
of loads, or the number of flights. The terms block and 
flight are used interchangeably. This entry is the number 
of times the list of loads is to be repeated. The next 
entry is the number of loads. The number of loads will 
always be the same as the quantity of cards 11 used as input 
The last entry is the number of types of loads. Loads are 
indicated as being a different type in order to alter the 
loading pattern during subsequent flights. It is important 
to understand that the program will always consider load 
types in ascending numerical order. A flight listing load 
types in the order 1, 2, 3, 4 will produce identical results 
with a flight listing load types in the order, 4, 3, 2, 1. 

8 . Card 10 

Only the limit load is listed here, with units of 
ksi. The maxima and minima of all cycles are input as a 
decimal fraction of this load entry. 
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9. Cards 11 



Information identifying each load cycle is presented 
on these cards. The first entry is an integer representing 
the step, or number of the load. This integer increments by 
one on each card, with the last card having a value equal to 
NLEVEL. Other entries on these cards are the type of load, 
the minimum and maximum decimal fraction of the limit load 
for the particular step, and the number of consecutive 
cycles at this load level. It is important to notice that 
the format requires the number of cycles to be entered as 
a real number, although fractions of cycles will not be 
processed . 

10. Cards NLEVEL + 11 

The first entry is an integer representing the number 
of the block, or flight. This entry is a dummy variable and 
is not used in any subsequent calculations. The other 
entries on this card indicate which types of loads are to 
be processed on this particular flight. Each type of load 
may be processed only once on each flight. The order in 
which the types are listed is of no significance, since the 
load types are automatically processed in ascending numerical 
order. 

B. RANGE -PAIR COUNTING METHOD 

In assessing the life consumed by an individual aircraft, 
damage is calculated from Miner's Law using S-N data. This 
type of life calculation is of greater value if the cycle 
counting method used takes into account the actual load-time 
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history of the aircraft. Many load cycle counting methods 
have been developed and used. Fatigue test experience indi- 
cates that a useful counting technique for aircraft fatigue 
must take into account the loading sequence as well as the 
magnitude of the maximum and minimum load peaks [Ref. 10] . 

The two most realistic counting techniques are the range- 
pair counting method and the rain-flow counting method. 

Both methods experimentally yield approximately the same 
results. The major' difference between the two methods is 
that the rain- flow method counts in terms of load ranges, 
or half cycles, whereas the range-pair method records the 
time history in terms of complete load cycles. When counting 
in terms of full cycles, the two methods are equivalent. 

The merit of counting half cycles is not important when 
analyzing aircraft fatigue because the large number of load 
reversals does not permit a half cycle to influence the 
numerical results. The range-pair counting method corresponds 
to the stable cyclic stress - strain behavior of a material 
in that strain ranges counted as cycles will form closed 
stress -strain hysterisis loops [Ref. 4], 

The computer program employed in this thesis has the 
option of using the range-pair counting method, or simply 
using each peak in the order that it occurs. When the range- 
pair counting method is not used, only momentary load values 
are considered, and information regarding the cyclic stress- 
strain pattern, which is important in fatigue calculations , 
is lost. Also, minor load variations, which are of little 
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importance in fatigue calculations, are counted as additional 
cycles . 

The range-pair counting method is illustrated in Figure 1. 
It counts a strain range as a cycle if it can be paired with 
a subsequent straining of equal magnitude in the opposite 
direction. For a complicated load history, some of the 
ranges counted as cycles will be simple ranges, such as 2-3, 
during which the strain does not change direction, but others, 
such as 1-8, will be interrupted by smaller ranges which will 
also be counted as cycles. In Figure 1 ranges are marked with 
solid lines and the paired ranges with dashed lines. 

Each peak is taken in order as the initial peak of a 
range, except that a peak is skipped if the part of the 
history immediately following it has already been paired with 
a previously counted range. If the initial peak of a range 
is a minimum, a cycle is counted between this minimum and 
the most positive maximum which occurs before the strain 
becomes more negative than the initial peak of the range. 

For example, in Figure 1 a cycle is counted between peak 1 
and peak 8, peak 8 being the most positive maximum that 
occurs before the strain becomes more negative than peak 1. 

If the initial peak of a range is a maximum, a cycle is 
counted between this maximum and the most negative minimum 
which occurs before the strain becomes more positive than 
the initial peak of the range. For example, in Figure 1 a 
cycle is counted between peak 2 and peak 3, peak 3 being 
the most negative minimum before the strain becomes more 
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Strain 




Figure 1 

Illustration of Range-Pair Counting Method 
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positive than peak 2. Each range that is counted is paired 
with the next straining of equal magnitude in the opposite 
direction, explaining why complete cycle rather than half 
cycle counts are made. For example, in Figure 1 part of the 
range between peaks 8 and 9 is paired with the range counted 
between peaks 1 and 8. 

C. DAMAGE CALCULATION 

Damage due to plastic strain and damage due to elastic 
strain is calculated separately. The behavior of the specimen 
is assumed to be elastic-perfectly plastic as illustrated in 
Figure 2. Behavior in compression is assumed to be the same 
as the behavior in tension. The total damage is determined 
using the Miner's summation D = Sn/N where 

D = total damage 

n = number of cycles at a particular load level 

N = number of cycles to failure at a particular 
load level 

Failure of the specimen is defined as occuring when D equals 
unity. 

Plastic strain damage is calculated using the Coffin- 
Manson theory [Ref. 5]. The Coff in-Manson theory is based 
on an experimentally determined log-log plot of plastic 
strain range vs. cycles to failure as shown in Figure 3. 

The slope of the line in Figure 3 is approximately -1/2 
for all metals, with the ordinate intercept making the 
behavior of each metal unique. The equation used for cal- 
culation of plastic strain damage is 
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stress , (S 



tensile yield 
stress, < 5 ty 






strain, e 

Figure 2. ELASTIC-PERFECTLY PLASTIC 
BEHAVIOR 





LOF 


log plastic 
strain range, 


\^slope approximately -$ 


6 P 






log cycles to failure, N 
Figure 3. COFFIN -MANSON DIAGRAM 
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N = ( e p / c ) k 

where N = number of cycles to failure at a specific plastic 
strain 

e p = magnitude of the plastic strain 

k = slope of the Cof fin-Manson curve, approximately 

- 1/2 

c = ordinate intercept; experimentally for each metal 
type 

For strain levels in the plastic range, failure occurs in 
4 

less than 10 cycles, and damage is calculated using the 
Cof fin-Manson theory. For strain levels in the elastic range, 
the Goodman diagram is the basis for calculation of damage. 

The maximum and minimum local strain levels are sequentially 
compared with second order least-squares curve fit data from 
the Goodman diagram. The curve fits used are for lives of 
10 4 , 10^, 10^, and 10^ cycles [Ref. 3], The elastic damage 
is totaled using the Miner summation and combined with the 
plastic damage, computed in a similar way, to yield the total 
damage. The amount of plastic damage and the total damage 
for a particular flight, and the total damage including all 
previous flights is printed at the end of each flight. 
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III. DISCUSSION OF RESULTS 



The load spectrum used in this analysis is the load spectrum 
A from MIL-HDBK-5A, and the loads are assumed to be applied at 
wing station 32 on an A7 aircraft wing [Ref. 2] . Load data 
were input as block loads and also as single randomized cycles. 
The Miner's Law summation of fatigue damage for this loading 
spectrum is .1300 per 1000 flight hours. Results obtained 
using this computer program are given in Tables II, III, and 
IV. It is important to notice that when fatigue damage is 
calculated without using the range-pair counting method, it 
makes only a minor difference on the total fatigue damage when 
the sequence of loading is changed or when negative loads are 
included . 

Block loads were arranged in LO-HI, HI-LO, and HI-LO-HI 
sequences. The results indicate the LO-HI sequence to be 
more damaging than the HI-LO sequence, with damage due to 
the HI-LO-HI sequence falling between the other two block 
load sequences. In the HI-LO sequence the HI loads leave 
a local residual compressive stress in the material, and 
the LO loads do less damage. 

In the case of RANDU generated load sequences, identical 
starter integers were used in order to evaluate the effect 
on the damage of inclusion of MIL spectrum A negative loads 
and ground cycles. The more negative loads that are included, 
the worse is the fatigue damage. The increase in damage with 
more negative loads is caused by the negative loads decreasing 
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TABLE II 



FATIGUE DAMAGE DUE TO POSITIVE MILSPEC A 
BLOCK LOADING, NORMALIZED TO 1000 HOURS 



SEQUENCE OF 
BLOCK LOADS 

LO-HI 

HI-LO 

HI - LO-HI 



USING RPCM* 
.16780102 
.069243014 
.10864216 

* RANGE -PAIR COUNTING METHOD 



NOT USING RPCM* 
.16780108 
.069243014 
.10839599 



TABLE III 

FATIGUE DAMAGE DUE TO MILSPEC A 



RANDOMIZED LOAD 


INPUTS, USING RPCM, 


NORMALIZED TO 


1000 HOURS 


INTEGERS USED 
TO START RANDU 


MINIMUM 


LOAD VALUES 






11 PERCENT 


NEGATIVE 


GROUND 




LIMIT LOAD 


MILSPEC A 


CYCLES 


83745,54711,54487 


.039226338 


.046267733 


.060848035 


13547,66549,7 


.095124505 


.085063167 


. 10429338 


9,583,4777 


.043321513 


.043640956 


.036781987 


48621,3,491 


.018672124 


.034077277 


.035881512 


73,559,1001 


.031819174 


.042006299 


.037130679 


357,833,1 


.037665060 


.049336702 


.050973855 
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TABLE IV 



FATIGUE DAMAGE DUE TO MILSPEC A RANDOMIZED LOAD 
INPUTS, NOT USING RPCM, NORMALIZED TO 1000 HOURS 

INTEGERS USED 



TO START RANDU 


MINIMUM 


LOAD VALUES 






11 PERCENT 
LIMIT LOAD 


NEGATIVE 
MILSPEC A 


GROUND 

CYCLES 


83745,54711,54487 


.040602274 


.047398917 


.049124956 


13547,66549,7 


.040137805 


.047727041 


.049612150 


9,583,4777 


.041204430 


.048346408 


.049540661 


48621,3,491 


. 041160211 


.049230121 


.050957575 


73,559,1001 


.040792227 


.047043338 


. 049311966 


357,833,1 


.041056350 


.049220584 


.050481893 
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the local residual compressive stress caused by positive loads 
into the plastic zone. This permits the local positve elastic 
stresses to be more damaging. 

Six plots are presented, beginning on page 33, to illus- 
trate the pattern of the individual load randomization process. 
Each plot portrays the first two flight hours of the aircraft. 
The range of values in the columns of Table III are caused by 
such variations in the load randomization pattern. 

The fatigue damage resulting from randomization of loads 
on a cycle basis, as shown in Table III, is significantly 
different from the damage due to the HI-LO-HI block loads 
illustrated in Table II. Vought Aeronautics Division con- 
ducted a fatigue study using randomized flights. Each flight 
was based on a mission profile specified by the Air Force. 
Fatigue damage resulting from the randomized flights was in 
close agreement with fatigue damage resulting from HI-LO-HI 
block loading. 
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IV. RECOMMENDATIONS FOR FUTURE RESEARCH 



Subsequent work in this area should be directed toward 
refining the input load sequence. The Naval Air Development 
Center collects accelerometer data from individual aircraft 
on a monthly basis. These load data could be used in this 
computer program to see if randomizing the load sequence over 
the life of the structure provides a different picture of the 
fatigue damage than randomizing by the month. Also, research 
is needed to determine the magnitude and sequence of the 
ground cycles. In this analysis all ground cycles were con- 
sidered of equal magnitude and equal spacing, which is not 
the case with a real airplane. A landing load spectrum should 
be used to be more accurate. 

The computer program itself could be improved by a better 
understanding and definition of the relaxation phenomena. 
Experiments to quantify the relaxation parameters will enhance 
the program. 

The Air Force is monitoring crack growth as the basis for 
their fatigue damage assessment. It would be productive to 
compare the total fatigue damage calculated by the two 
methods . 
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onnoonooooonoonoonooooooooooooooonooonnnoooooooonooooooooooooooonooo 



COMPUTER PROGRAM 



MODULE I 

INPUT ROUTINE FOR THE 
SEQUENCE ACCOUNTABLE FATIGUE ANALYSIS 



INPUT 

DATA CARD 1. NDECK = THE NUMBER OF DATA CECKS TC 

BE RUN SEQUENTIALLY 
I PR I NT= THE VALUE CONTROLLING THE 
WRITE STATEMENTS 

1. PERMITS MAXIMUM PRINTOUT 

2. SUPPRESSES RANGE-PAIR 
PRINTING 

3. MAXIMUM SUPPRESSION CF 
PRINTOUT 

I R PCM = ThE VALUE CONTROLLING THE 
ENTRY INTO THE RANGE-PAIR 
COUNTING SUBROUTINE 

1. ENTER RANGE-PAIR COUNTING 
SUBROUTINE 

2. SKIP RANGE-PAIR COUNTING 
SUBROUTINE 

FORMAT 3 I A 

EACH CATA CECK CONTAINS THE FOLLOWING CARCS - 
CARD L. 



CAPC 2. 



TEST IDENTIFYING INFORMATION 
FGRMAT 1 6 AA 

TM = MATERIAL TYPE 
TYS = TENSILE YIELD STRESS (KSI) 
EPSD = L CF STRAIN INTERCEPT 
CQF MAN =1 NVERSE CF COF F I N-MANSQN 
SLOPE 

ELMCD = MODULUS CF ELASTICITY (KSI) 
FORMAT AAA,3F18.5, FIO. 2 



CARCS 3 , ,6. 



A 


(N) 


N=A 


» 7 ( A,e,C AR 


E C 


L - 


FF 


IC 


I ENT 








OF SECOND 


CR 


CE 


R 


LE 


AST 


B 


(N) 




SQUARE FI 


T 0 


F 


5- 


N 


CATA 








FOP CURVE 


CF 


I 


0# 


*N 




C 


(N) 




CYCLES.) 


















(STMAX = 


A ( N 


)* 


STM I 


N**2 








+ B(N ) 


*ST 


MI 


N 


+ 


C ( N ) 


T 


ITL 


El 


(TITLElt T 


ITL 


E2 


* TI T 


LE3, 


T 


ITL 


E2 


TITLEA ID 


ENT 


IF 


IE 


S 


THE 


T 


ITLE3 


SOURCE OF 


TH 


E 








T 


ITL 


EA 


S-N DATA) 












N 






(PUNCHED 


IN 


CC 


LU 


MN 


72. 








FOR INFO. 


) 










MATE 


RIAL 


TYPE (COLUMNS 


73 


• 


80 


• 


FOR 



I NFG ONLY. ) 



CARC 7. 



FORMAT 


3F i 8 . 5 , AAA 


Cl 


(CONSTANTS TO BE USED I 
CALCULATION OF ECLILIBR 


El 


PERIOD, ENEP.) 

( ENEP=C1/(KTSMAX**E1*KT 


E2 


**E2 ) ) 


FORMAT 


3F18.5 
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CARD 


8. 


AKT = 


STRESS CONCENTRATION FACTOR 
USED THE FIRST TIME THROUGH 








THE PROGRAM 






FORMAT 


F18.5 


CARD 


9. 


NBLOCK = 


NUMBER OF BLOCKS (NO. OF 
TIMES TO REPEAT LIST OF 
LOADS ) 






NLEVEL = 


NUMBER OF L3ACS 






NTYPE = 


NUMBER OF TYPES CF LOADS 






FGRMAT 


3 110 


CARD 


10 . 


TLL = LIMIT LOAD (KSI) 






FORMAT 


F18.5 


CARDS 


1 1 » • * • t 


NLEVEL + 


10. 






K 


THE KTH STEP (K=l, NLEVEL) 






ITYPE (K) 


= IDENTIFYING TYPE OF 








THE KTH LOAD 






STMIN(K) 


= THE KTH MINIMUM (DECIMAL 








FRACTION OF TLL) 






STMAX (K) 


= THE KTH MAXIMUM (DECIMAL 








FRACTION OF TLL) 






ENN (K) 


= NUMBER OF CYCLES AT THE 








KTH LOAD 






FORMAT 


2(I4,2X),3(F18.5f IX) 



CAROS N LEVEL + 11 ,..., NLEVEL + 10 + NBLGCK . 

JJ THE JJTH BLOCK ( J J = l , NB L CC K ) 

NN ( J J » 1 ) = TYPE CF LOAC INCLUDED IN 
JJTH BLOCK 

NN(JJ,2) = (THERE WILL BE CNE NN 
VALUE ON THE C ARC FOR 
NN(JJ, ) = EACH DIFFERENT TYPE CF 

LOAC INCLUDEC IN THE JJTH 
BLC C K ) 

FORMAT 915 
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CCMMCN/MSAL/RNCYC ( 42 05) ,KPMAX, I PR INT 

C PENSION RTMIN ( 4205) »RTMAX(420 5) , ITYPE( 420 5) ,RNN'(42C5 

L) 

CCMMGN/MCCR1/NLEVEL, IRPCM, I PLOT t E LMOD , T Y S , E PS C , CO FRAN , 
1RES( 4205) ,AKT,SUMENN, SUMNC,C1, E2,E1 
CCMMCN/MCCR2/STMIN( 4205),STMAX( 4205), ENN( 4205), 

1C( 10) ,R( 10),NN<20, 10) ,A(10) ,B( 10) 

C IRENS IGN HOLD(lOl) 

C 

C 

C INSERT RANDU PACKET HERE - REMOVE ISKIP=0 CARO 

I S K I P = 0 

P E AC ( 5 , 4 ) NDECK, IPRINT, IRPCM 

4 FCFMAT (314) 

WF ITE (6,3) NDECK 

3 FCPMAT( 1H1 , I4,32H DATA DECKS ARE TO BE PROCESSED.) 

IF (IRPCM. GE. 2) GC TO 6 
WFITE(6,U) 

GC TO 13 

6 V F I TE ( 6 , 7 ) 

7 FORMAT (25H NO COUNTING METHODS USED) 

13 CONTINUE 

CC 595 LMN = 1, NDECK 
Vi F ITE ( 6 ,5 ) 

5 FORMAT (40H SEQUENCE ACCOUNTABLE FATIGUE EVALUATION) 
FEAD(5,5)T1,T2,T3,T4,T5,T6,T7,T8,T9,T10,T11 ,712,T13,T1 

14 , T 15 , T 16 

9 FCFMAT ( 16A4) 

WPITE(6,8)T1,T2, T3,T4,T5,T6,T7,T8,T9,T10 ,T1 1,T12,T13,T 
114 ,T15,T16 

8 FCFMAT (16H1SPECTRUM FROM ,16A4) 

C 

C 

C INPUT OF DATA PECULIAR TC A MATERIAL 

C 

C 

P E AC ( 5 , 10)TM1,TM2,TM3,TM4,TYS,EPSD,C0FMAN,ELMCD 

10 FCFMAT (4A4,3F18.5»F10.2) 

F I TE ( 6» , 12)TM1,TM2,TM3,TM4,TYS, EPSCtCOFMAN, ELMOD 
12 FC FMAT ( 18H MATERIAL TYPE — ,4A4,//16H TENSILE YIELC S 
1141-TReSS (KSI) — ,F18.5, 

1 //23H LCF STRAIN INTERCEPT =,F18.5,// 

' 13 1H INVERSE OF C OFF I N-MANSON SLOPE ,F 18 . 5 , // 

1 1 £F ELASTIC MODULUS =,F18.5) 

PEAD(5, 14) (A(N) ,3(N) ,C(N) ,TITLE1,TITLE2, TITLES, TITLE4, 
IN =4 ,7 ) 

14 FCFMAT (3F18.5 ,4 A4 ) 
kFITE (6,16) 

16 FCFMAT (44HOCOEFFICIENTS OF SECOND ORDER LEAST SQUARE F 
114HIT OF S-N DATA) 

WRITE (6, 18) 

L8 FCFMAT (43H SMAX = A ( I ) *SM IN**2 + B(I)*SMIN + C(I)) 

WFITE(6,20) 

20 FORMAT ( 5X ,5H LIFE,10X,5H A(I),14X,5H 3(I),14X,5H C(I)) 
WFITE(6,22) (N,A(N) ,B(N) ,C(N) ,N=4,7 ) 

22 FORMAT { SH 10** , 1 2 , 3F 18. 5 ) 

WRITE(6,56)TITLE1,TITLE2,TITLE3,TITLE4 
56 FORMAT ( 39H0UNN0TCHED COUPON S-N CATA CERIVEC FROM/ 

12SH INFORMATION SUPPLIED FROM ,4A4) 

REAC(5,14) C 1 , El , E 2 
VF ITE (6 ,24) 

24 FORMAT (36H OR ESI DUAL STRESS RELAXATION FUNCTION) 

VFITE(6,26) 

26 FCFMAT ( /41H ENEP = C 1 / ( KTSM A X**E 1 * KT S M E AN** E 2 ) / ) 

WFITE(6,28) C 1 , E 1 , E2 

28 FCFMAT ( 13H+WHERE Cl =,E15.8,9H , El =,Fi0.3, 

1 1CF ANO E2 = , F 1 0 • 3 ) 

C 

c 
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INPUT OF DATA PECULIAR TC A SEQUENCE 



C 
C 
C 
C 

R E AC ( 5 » 65 ) AKT 
65 FORMAT! F18.5 ) 

R£AC(5,32) NBLOCK, JLEVEL, NTYPE 

32 FORMAT (3110) 

l« F ITE (6 ,34 ) NBLOCK , J LEVEL 

34 FCFMAT (//I 10, 23H TIMES THROUGH BLCCK OF, 110,61- LOADS) 

R EAC( 5» 35 ) TLL 

35 FCFMAT! F18.5) 

VsP ITE ( 6 ,33 ) TLL 

33 FCFMAT ( / 5X,13H LOAD LIMIT =,F18.5) 

IF! ISKIP.EG.100) GC TC 37 

R EAC(5, 83 ) ( IOUMMY , ITYPE(K) ,RTM IM K ) ,R T M AX ( K ) , RN'N ( K ) , 
1K=1 , JLEVEL ) 

88 FCFMAT (I4,2X,I4,2X,F18.5,1X,F18.5,1X,F18.5,1X) 

37 hFITE!6,38) 

38 FCFMAT ( // i 1H STEP TYPE,10X,6H STMIN,14X,6H STMAX,15X, 
141- 6NN ) 

IFdSKIP.EC. 100) GG TO 89 

hFITE(6»36) (K,ITYPE(K), RTMIN(K) , RTM AX ( K ) ,RN N ( K ) , K= 1 , 
1JLEVEL) 

GC TC 91 

89 fcFITE!6,36)IK,ITYPEIK) , R T M I N ( K ) ,RTMAX(K) ,RNN(K),K = 2C1, 
1JIEVEL, 250 

91 b F ITE <6,39) 

36 FCFMAT! IX, 14, 2X, I4,2X,F18.5,1X,F18.5,1X,F18.5) 

39 FORMAT ( //47H BLOCK TYPE TYPE TYPE TYPE TYPE TYPE TYPE 

15F TYPE/) 

CC 42 JJ=1, NBLOCK 

R EAC ( 5 , 40 ) IDUMMY , ( NN ( J J , KK ) , KK= 1 , NT Y PE ) 

V> R I T E < 6 ,40 ) JJ , INN ( J J » KK ) ,KK=1,NTYPE ) 

40 FCFMAT! 1315) 

42 CC NT I NU E 

SO MENN= 0. 0 

SLMNC=0.0 

PES!1)=0.0 

WRITE!6,8)T1,T2,T3,T4,T5,T6,T7,T8,T9,T10,T11,T12,T13,T 
114, T15, T16 
UFITE(6,51) ( AKT ) 

51 FC FMAT ( // 7H AKT = ,(F6.2)) 

W F ITE ( 6 » 55 ) Cl 

55 FCFMAT ( //24H RELAXATION CONSTANT C1=,F15.2) 

IF ! IR P CM . GE . 2 ) GC TO 59 
to F ITE ( 6 , 11 ) 

11 FORMAT (38H0SPECTRUM SUBJECTED TC THE RANGE-PAIR 

118FCCUNT ING TECHNIQUE) 

59 CONTINUE 

DC 1002 KFL=1, NBLOCK 
JJJ = 1 

CC 60 J=l, JLEVEL 
CC 70 K K= 1 , NTYPE ‘ 

IF !NN( KFL,KK) .EQ.O) GO TO 60 
IF ( ITYPE(J) .EQ.NN!KFL,KK) ) GO TC 150 
70 CONTINUE 

150 STMIN! JJJ)=RTMIN( J)*TLL 
STM AX! JJJ )=RTMAX( J)*TLL 
ENN! JJJ )=RNNi J ) 

J V J = J J J +1 

60 CONTINUE 

N LEVEL = J J J-l 
CALL CORE ! KFL ) 

IFR INT = 2 
1002 CONTINUE 
597 CONTINUE 

595 CONTINUE 

596 CONTINUE 
580 STOP 

ENC 

C 
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SLERGUTINE CORE ( KFL ) 

C 

C CORE PROGRAM OF 

C THE SEQUENCE ACCOUNTABLE FATIGUE ANALYSIS 

C 

c 

c MODULE II 

C 

C LOCAL STRESS AND STRAIN DETERMINATION 

C 

C 

CCMNCN/MCEC1/SIGMAX( 4205 ) , S IGM IN ( 4205) 

COMMON/ MS AL/RNCYC ( 42 05 ) t KPMAX , I PRINT 
DIMENSION PLSTRA ( 4205), EN( 4205), EX( 4205) 

COM MC N/ MCORi/NLEV EL, I RPCM, I PLOT , ELMOD , TY S , E PS C , CC FMAN , 
IRES ( 4205) ,AKT,SUMENN» SUMNC ,CI,E2,E1 
CCMMGN/MC0R2/STM INI 4205),STMAX( 4205), ENN( 4205), 

IC (1C) ,R(10) ,NN< 20,10) ,A(10) ,B(10) 

J 2 = K F L 

IF (JJ.GT.3) I PR I NT = 3 
I R A IN=1 

WRITE(6,54) JJ 

54 F C R MAT ( 20H FLIGHT OR BLOCK NO., 15) 

IF ( IPRINT.GE.3) C-C TO 61 
to F 1TE(6,62 ) 

62 FC RMAT ( IX ,6H STMAX , 3 X , 5 HSTM IN ,3 X ,6FS I GMA X, 3 X , 6FS I GM I N , 

13X,2HRES,3X,5HEQRES,3X,3HENN,11X,3HNEP) 

61 CONTINUE 

C 

C 

148 OC 570 J=1,NLEVEL 

149 I=J+1 
PLSTRA ( J ) =0. 0 

C 

c 

C DETERMINE SEQUENCE ACCOUNTABLE RESIDUAL STRESS AND 

C 

C CORRESPONDING PLASTIC STRAIN 

C 

150 ASMAX=AKT=*$TMAX( J) 

ASMIN=AKT=*STMIN(J) 

ASM EAN = ( ASMAX+ASM IN ) / 2 • 

JA= C 

170 IF (RES ( I- 1 ) +ASM I N + TY S ) 190,190,200 
190 JA= 1 

AAA=RES ( I-D+ASMIN 

PLSTRA(J) = -L.=MAAA/ELMOD)*( l.+AAA/TYS) 

200 JE = C 

202 IF (RES(I-1)+ASMAX-TYS) 220,210,210 

210 J E=- 1 

A AA =A SMAX 
EEE=TYS-RES( I — 1 ) 

IF (BBB.GE.TYS) GO TO 214 
212 E E E = 0. 0 

214 FLSTRA(J) = AAA*AAA/ ( ELMOD«TYS )-EEB*BB9/ ( ELMOC=»TYS ) 

220 IF (JA+JB) 230,250,240 
230 RES(I)=TYS-ASMAX 
GC TO 290 

240 RES( I )=-TYS-ASMIN 
GC TO 290 

250 IF(JA) 260,260,270 
260 FES(I)=RES(I-1) 

GC TO 290 

270 R E S ( I ) =-ASMEAN 

290 S IGMAX( IRAIN ) =RES ( I ) +A SMAX 
S I GM IN ( IRAIN) =RES ( I J + ASMIN 
RNCYC ( IRAIN) =ENN ( J ) 
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IF ( AS MAX . LE . TYS ) GO TO 410 
400 ECRES=-ASMAX+TYS 
GC TO 440 

410 IF (ASMIN.GE.-TYS) GO TO 430 
420 ECRES=-ASMIN-TYS 
GC TO 440 
430 E OR ES = 0 .0 

440 CIF=RES(I J-EQRES 
C 
c 

C CALCULATE RELAXATION FUNCTION 

C 

C 

ABNAX=ABS(ASMAX) 

AEMIN= ABS (ASMIN ) 

AENEAN=ABS(ASMEAN) 

IF (ABMAX.LT. I. ) ABMAX =1 . 

IF (ABMEAN.LT.l.) ABME AN=0. 5 
IF ( ABM I N-ABMAX ) 444.444,442 
442 A£M=ABMIN 
GC TO 446 
444 A EN=ABMAX 

446 ENEP=Cl/(ABM=*~El : i e A9MEAN~ 5 *E2) 

IF (IPRINT.GE.3) GG TO 351 

WPITE(6,350) STMAX(J) tSTMIN(J) .SIGMA X(IRAIN) ,SIGMIN(IR 
IAIN),RES( J), EQRES,ENN(J), ENEP,J, I RAIN 

350 FORMAT ( 6 ( F7. 2 ,1 X ) , F6. 2 , IX , E 15 . 8 , 1 5 , 1 5 ) 

351 CONTINUE 
C 

C 

C CALCULATE RESIDUAL STRESS RELAXATION 

C 

C 

IP AIN=I RAIN+1 
A EC I F = A BS ( Cl F ) 

C-C TO 360 
C 

C CMITTED "GO TO 560" CARD FROM THIS POSITION, AS ShC'*N 

C IN MANUAL 

C 

360 IFCABOIF.LT. 5.) GC TO 560 
370 IF (1000.*ENN< J) .LT.ENEP) GO TO 56C 
IF ( ENN (J).LE.IO.) GO TO 560 
380 NFLAG=0 

N F L AG2= 10 
IRAIN=I RAIN-1 
CLMMY=ENN( J ) 

450 IF ( DUMMY-ENE P ) 470,460,460 
460 C LMMY = DUMMY/2 . 

NF LAG=N FL AG + 1 
GC TO 450 

470 CYCINT=CUMMY/10. 

CC 500 K= 1 , 10 
C ECK=FLGAT ( K ) 

EN (K)=CYCINT*DECK 
IF (K.EC.l) GO TO 490 

480 EX(K)=EXP(-2.303*EN(K-1)/ENEP)+EXP<-2.303*EN(K)/ENEP) 

GC TO 500 

490 E > (K ) =1 . + EXP (— 2 .303*EN ( K ) / ENEP ) 

500 CONTINUE 

IF(NFLAG.EC.O) GO TO 530 
510 NFLAG2=NFLAG+10 

CC 520 K=ll , NFLAG2 
E N ( K ) =2 . ^DUMMY 

EX(K)=EXP(-2.303*EN(K-1)/ENEP)+EXF(-2.303*EN(K)/ENEP) 
CU M MY = 2 . ^ DUMMY 
520 CONTINUE 
530 CC 559 K=1 ,NFLAG2 

S IG'MAX ( IRAIN)=ASMAX+ECRES+DIF*EX(K)/2. 

S I GM IN ( IRAIN)=SIGMAX( I RAIN ) -ASMAX+ASMIN 
RNCYC( I PA IN ) = EN ( K ) 

545 IF (K.EG.l) GO TO 540 
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539 RNCYC { IRA IN ) = RNC YC ( IRA I NJ -£ N< K- 1 ) 

540 IF (IPRINT.GE.3) GO TO 55i 
WRITE(6»550)SIGMAX(IRAIN) , S IGMIN ( IRA IN) t RNCYC ( I RA IN ) , I 

1RAIN 

550 FCPMAT( 16H RELAXATION , 2 (F7 . 2 , IX ) , 1 6X , F6 .2 ,3 1 X , I 6 ) 

551 CONTINUE 

I P A 1N= I R A I N+ 1 

559 CONTINUE 

560 CONTINUE 

RES (I )=EQRES+DIF*EXP(-2.303*ENN( JJ/ENEP) 

570 CONTINUE 
569 RE$(1)=RES(J) 

IN=IRAIN-1 

C 

C 

C 

C MODULE III 

C CYCLE COUNTING TECHNIQUE 

C 

(3 45$$#%****%*#*****$***$*###*#*# If ****%*&■$** $■$****$****# 

C 

c 

C CALL SUBROUTINE TO RANGE-PAIR COUNT SPECTRUM 

C 

C 

I F ( IRPCM.GT .1 ) GO TO 591 
CALL RPCM ( IN) 

GC TO 592 

591 CCNTINUE 
KFMAX=I N 

592 CCNTINUE 
C 

C 

(3 * * * *#### ###* #**####* :* * :* $£ *:£ ######## ##zf ######**#*#*### # 

C 

C MODULE IV 

C DAMAGE ACCUMULATION CALCULATION 

C 

(3 * ***$$# **$****##$ ****** ******** *#*##**#** ** *** ******** 

C 

C 

IF ( IPRINT.GE.3) GO TO 552 
WFITEC6 ,53) 

53 FORMAT (//1X,39H LOCAL STRESSES ANC PLASTIC STRAINS WIT 
124HH RESULTING FATIGUE L I F E//10X , 4HSTEP , 10X , 
114FPLASTIC STRAIN, 10X,10HMAX OR M IN , 15X , 6F0 AMAGE ) 

552 CCNTINUE 
C 

C 

C CALCULATE DAMAGE FROM PLASTIC STRAIN CYCLES 

C 

C 

SUM0EL = 0 . 

DC 531 JKL=1,NLEVEL 
A A = 1 . 

IF (PLSTRAt JKL) ) 532,531,533 

532 A A = - 1 . 

533 PLSTRA( JKL)=AA*PLSTRA(JKL) 

CYCLES= ( PLSTRA( JKL )/EPSD) **CCFMAN 
C AMI. /CYCLES 
SLMNC=SUMNC+OAM 
SLMCEL=SUMDEL+DAM 

I F ( I PR INT . GE .3 ) GO TO 531 
IF (AA) 535,535,537 
535 WRITE (6, 199) JKL, PLSTRAt JKL) ,CAM 
199 FCFMAT ( 10X ,14 ,1 2X , F 10 . 5 , 15X , 3HM I N , 10X , E 14 .6 ) 

C-0 TO 5 31 

537 WFITE(6,219) JKL,PLSTRA(JKL) , CAM 

219 FCPMAT(10X,I4,12X,F10.5,15X,3HMAX,10X,E14.6 ) 

531 CCNTINUE 

WFITE(6,541) SUMDEL 
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541 

13 

536 

C 

C 

C 

C 

C 

C 

534 

310 

447 

320 

330 

332 

334 

335 

336 

338 

339 

340 

599 

600 

593 

575 

C 



FORMAT ( 35X » 29H DAMAGE FROM PLASTIC STRA I NS= , E 1 5 . 8 ) 
IF ( IPRINT.GE.3) GO TO 536 
WR ITE (6 , 13 ) 

FORMAT ( / 8X,15H SIGMAX SIGMIN, 8X,6H RNCYC * 8X, 
12SF CYCLES ENN/CYC) 

CONTINUE 



CALCULATE ELASTIC CYCLE DAMAGE FROM LEAST SQUARE 

FITTED S-N DATA (MODIFIED GOODMAN DIAGRAM FORMAT) 



DC 600 JKL = 1 * KPMAX 
T 1 Y S=TY S/ 5 . 

IF(SIGMAX(JKL)-SIGMIN(JKL).LT.1.6*TYS) GO TC 310 
CYCLES=10.**4. 

GO TO 340 

1F(SIGMAX( JKD.GE.TTYS) GO TO 320 
CYCLES= 10 .**9. 

C-C TC 340 
N =4 

CC 330 M= 1*4 

R (N)=A( N)*SIGMIN( JKL )**2 + B ( N ) *S I C-M IN ( JKL )+C ( N )-S I GM a X ( 
1JKL) 

N = N + 1 
CONTINUE 

IF(R(7)*R(4) ) 338,338,334 
AER7=A8S(R(7) ) 

A E R 4=A 8 S ( R ( 4 ) ) 

IF(A8R7.LE.ABR4) GO TO 336 
EXFC=4.+R(4)/(R(4)~R(5) ) 

GO TO 339 

EXFC=7.+R(7)/(R(6)-R(7) ) 

C-C TO 339 

S LMR=R ( 4)+R( 5)+R(6 )+R(7 ) 

SLMR2 = P (4)**2+R< 5 )**2+R(6 )**2+R ( 7)**2 
SLMR3=R(4)**3+R( 5 ) **3 + R ( 6 ) **3+R (7)**3 
SLMR4=R (4)**4+R( 5 )**4+R (6)**4+R ( 7)**4 
SUM FN=4 .^R (4}+5.*R(5 )+6.*R(6)+7.*R(7) 
SLMR2N=4.*R(4)**2+5.*R< 5)**2+6.*R (6)**2+7.*R ( 7) **2 
CE L 1=4. *SUMR2*SUMR4-4. *SUMR 3**2 
CEL2=SUMR*SUMR2*SUMR3-SUMR4*SUMR**2 
CEL3=SUMR*SUMR2*SUMR3-SUMR2**3 
CC1=22.*SUMR2*SUMR4-22.*SUMR3**2 
DC2=SUMR2*SUMR3*SUMRN-SUMR*SUMR4*SUMRN 
DC2=SUMR*SUMR3*SUMR2N-$UMR2N*SUMR2**2 
EXFC=(DC1+DG2+D03 )/ ( DEL 1 + DEL2+DEL3 ) 

CYCLES=10.**EXP0 
IF(EXP0.LE.4. )CYCLE$=10.**4. 

ENNCYC=FNCYC (JKL) /CYCLES 
SUMNC=SUMNC+ENNCYC 
SLMCEL= SUMDEL+ENNCYC 
I F ( IPR I NT . GE .3 ) GO TG 600 

WRITE (6, 59 9) SIGMAX (JKL) , SIGMIN (JKL) ,RNCYC( JKL) , CYCLES 
1 , E N NCYC 

FO FM AT ( 3X,2(F7.2,1X) , 6X,F6.0, 7X ,2 ( El 5 . 8 , IX ) ) 
CONTINUE 

WFITE (6,593) SUMDEL 

FORMAT ( /35X,21H DAMAGE PER THIS SET=,E15.8) 
WRITE(6,575) SUMNC 

FOFMAT ( /35X , 18H TOTAL ENN/CYC =,E15.8) 

RETURN 

ENC 
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SUBROUTINE RPCM(NPKS) 

Tt-IS PROGRAM EMPLOYS THE RANGE-PAIR CYCLE 
METHOD TO GENERATE ANALYSIS SPECTRUM FROM 
LC AC SPECTRUM 



COUNTING 
A GIVEN 



PROGRAM ARRAYS 
( INFORMATION NEEDED TO 



CHANGE DIMENSIONS) 



ARRAY NAME 
SIGMAX 
KK 



SIGMIN 

FNCYC 

NSTEP 

RES 

INDEX 

CYCLE 

RNECYC 

NNSTEP 

ISAVc 



DEF INITION 

PEAKS OF THE INPUT LOAD SPECTRUM 
THE NUMBER OF ADDITIONAL CYCLES 
(EXCLUDING INPUT CYCLES) WHICH 
THE PROGRAM WILL GENERATE 
VALLEYS OF THE INPUT LOAD 
SPECTRUM 

K COUNTERS OF THE PEAKS AND 
VALLEYS 

STEP NUMBERS OF THE INPUT 
SPECTRUM 

RESIDUE SPECTRUM 
STEP NUMBERS OF ELEMENTS IN RES 
RANGE-PAIR CGUNTED CYCLES 
K COUNTERS OF THE CYCLES OF THE 
UNSORTED ANALYSIS SPECTRUM 
STEP NUMBERS OF THE ELEMENTS CF 
THE UNSORTED ANALYSIS SPECTRUM 
VALUES OF NSTEP(J) SUCH THAT 
RNCYC(J) IS <1.0 ANC VALUES OF 
NSTEP(J) SUCH THAT SIGMAX(J-I)= 
SIGMAX(J) ANC SIGMINl J-1)=SIGMIN 
CC MMCN/ MS AL/ RNCYC ( 4205 ) , KPMAX , I PR I NT 
CCMMCN/MDEC1/SIGMAX( 4205 ) , S I GM I N ( 4205) 
CCMM0N/MCEC2/NSTEP ( 4205) , LR , KMA X , KM IN , K 3 1 
COMMGN/MDECR/RES ( 4205),INDEX( 4205)»INDltI 
104, KIND 

CCMMQN/MCYG/CYCLE(200,2),RNECYC( 4205) , NNST 

COM MON/MCGDE/L, LIND 

DIMENSION ISAVE(4205) , TITLE (8) 

9959 N PUNCH = 0 

CC 6000 1 = 1, N PKS 
80CC N S I EP ( I ) = I 

IF ( IPRINT.GE.2JG0 TO 103 
WRITE(6,20)NPKS 

2C FCRMAT( 1H0,40HTHE NUMBER OF PEAKS OR VALLEY 
1 2CH PUT LOAD SPECTRUM = ,15//) 

WP ITE(6 ,22) 

22 FORMAT (3SX,5HSIGMA/10X,4HSTEP, 1 3X , 7HMAX I MUM 
17HMINIMUM,12X,9HCCUNTER K/ ) 

WP ITE(6,25) ( NSTEP ( I ) , SIGMAX ( I) ,SIGMIN( I ) »R 
l 1 , NPKS ) 

25 FC PMAT ( 8X , 



DIMENSION 
NPKS + KK 



NPKS + 

NPKS + 

NPKS + 

2*NPKS 
2*NPKS 
NPKS + 
NPKS + 

NPKS + 
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KK 

KK 

KK 



KK 

KK 

KK 



( J ) 

ND2, IND3, IN 
E P ( 4205 ) 



S IN THE IN 



03 



CONTINUE 



,16X, 

NCYC ( I ) , I = 
15 ,10X,E13.6,10X,E13.6,10X,F10.5 ) 



SORT THROUGH THE LOAD SPECTRUM - FULL OUT THOSE PEAKS 
ANC VALLEYS WHERE COUNTER K IS LESS THAN l.C 



J = 1 
L = C 

NFES = 1 
NCYNO = 100 
J M A X = 0 
CC 100 1=1, NPKS 

IF (RNCYC ( I ) .GE. 1.0) GO TO 100 
XI = SIGMAX(I) 

X 2 = SIGMINM I ) 

CALL CYCGEN(X1,X2,RNC YC( I) ,NSTEP( I) ) 
ISAVE(J) = I 
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j = j +1 

100 CCMINUE 
JMAX = J - 1 

NFKSN = NPKS - JMAX 
IF (JMAX. EG. 0) GO TO 200 
WFITE(6,23) (ISAVE(K),K = 1,JMAX) 

IF( IPRINT.GE.2JG0 TO 101 

23 FCPMAT( 1H0,40HSTEP NUMBERS OF THOSE PEAKS A NC VALLEYS 
1 5 IF IN THE LOAD SPECTRUM WHOSE COUNTER K IS LESS THAN 1 
12F.0//( 1717) ) 

101 CONTINUE 

DC 110 J = 1 » JMA X 
I = ISAVE(J) - (J-l) 

NFKN = NPKS - J 
I F ( I.EQ.NPKN) GO TO 110 
CC 115 II = I, NPKN 
S IGMAX( II ) = S IGMAX ( 1 1 + 1 > 

SIGMINi II) = S I GM IN ( 1 1 + 1 ) 

N STEP C II) * NSTEPUI + 1) 

PNCYC(II) = RNCYC ( 1 1 + 1 > 

115 CONTINUE 
110 CONTINUE 
2CC CONTINUE 



SCPT THROUGH THE LOAD SPECTRUM DATA-C0M3 IN E STEPS WITH 
IDENTICAL PEAKS AND VALLEYS WHICH CCCUR CONSECUTIVELY 



J * 1 

DC 300 I = 2 » NPK SN 

I F ( S IGMAX ( I ) .NE. S IGMAX C I— 1J > C-C TO 300 
IF (SIGM IN( I ) .NE. SIGMINt I — 1 ) ) GC TO 300 
ISAVE(J) = I 

PNCYC(I-l) = RNCYC ( I- 1 ) + RNCYC ( I) 

J = J + 1 
300 CONTINUE 

IF (J.EC.l) GO TO 6000 
JMAS = J - 1 
CC 311 J = 1 » J M AS 
I = ISAVE(J) - (J-l) 

NFKN = N'PKSN - J 

IF (I . EQ. NPKN) GO TO 311 

CC 316 II = I, NPKN 

S I G M A X ( II ) = S IGM AX ( 1 1 + 1 ) 

SIGMINt II) = SIGMIN(II+i) 

NSTEP(I I ) = NSTE P ( II + l) 

PNCYC(II) = RNCYC ( 1 1 + 1 ) 

316 CONTINUE 
1 • CONTINUE 

NFKSN = NPKSN - JMAS 



RANGE PAIR COUNT THE ADJUSTED LCAC SPECTRUM 



600C 1=1 

Ke = o 

L= JMAX 
KM IN=0 
KM AX=0 
LP =0 
K 3 1 = 0 

1 IF (RNCYC(I) .GT. 1.0) GO TO 400 

IF (KB .NE. 0) GO TO 5 
X 1=S I GM AX ( I ) 

X2 = S IGM IN ( I ) 

IND1 = NSTEP( I ) 

INC2=IND1 
1 = 1 + 1 
KE=1 
GC TO 1 
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5 >3=SIGMAX( I ) 

X4=SIGM IN ( I ) 

I NC2 = NSTE P ( I ) 

INC4=IND3 
KM IN=1 
KMAX=0 
K2 1=0 

IF (RNCYC(I) .EQ. 1.0) C-0 TO 6 
KEY = 1 
K I N D= 1 
GC TO 415 

6 K E Y=0 

CYCNC=A INT(RNCYC( I ) +0.5 ) 

CALL CECIDE(X1 ,X2 ,X3,X4 t KEY, I t CY CNO » KCY GEN ) 

100C GC TO ( 10,10,30) ,KCYGEN 

10 KE=1 

C 
C 

C REVERSE ORDER OF NEXT TWC CARDS TC RUN. "1=1+1" GCES 

C AFTER " I F< KMIN.NE . 1) GO TO 36" AS PER ERRATA SHEET IN 

C PRINTED MANUAL. CHECK PROGRAM WITH CARDS IN THE 

C MANUAL'S CRDER 

C 
C 

IF (KMIN .NE. 1) GO TO 36 
1 = 1 + 1 

IF (I .LE. NPKSNJ GO TO 5 

RES (LR+ 1) = XI 

RES (LR+2 ) = X2 

INCEX( LR+1 ) = I ND1 

INDEX ( L R + 2 ) = I ND2 

LR MAX = LR+2 

GC TO 2000 

30 IF (KMIN .NE. 1) GO TO 35 

12 1=1+1 

IF (I .LE. NPKSN) GO TO 31 
R ES ( LR+ 1 ) = XI 
RES (LR+2 ) = X2 
R E S ( LR+2 ) = X3 
I NCEX ( L R+ 1 ) = IND1 
INCEX( LR+2 ) = I ND2 
INCEX ( LR+3 ) = I ND3 
LRMAX = LR+3 
GC TO 2000 

31 X4 = S I GM A X ( I ) 

IND4=NSTEP( I ) 

KM A X= 1 

KM I N=0 
K2 1 = 1 

32 IF (RNC YC ( I ) .GT. 1.0) GO TC 40 
GC TO 6 

40 KEY = 1 

KINC = 0 
GC TO 415 

35 X4 = SI GM IN ( I > 

I NC4 = NSTEP(I) 

KMIN = 1 
K M AX = C 
K21 = 0 
GC TO 22 

36 X3 = SI GMIN( I ) 

INC3 = NSTEP(I) 

KMIN = 1 

KMAX = 0 
GC TO 12 
400 KEY = 1 

IF (KB .NE. 0) GO TO 410 
XI = SIC-MAX (I) 

X2 = S I GM IN( I ) 

X 2 = SI GMAX( I ) 

X4 = SI GM I N( I ) 

INC1 = N STEP ( I ) 
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I NC2 = IND1 
INC3 = IND1 
INC4 = IND1 
KMIN = 1 
KNAX = 0 
K2 1 = 0 

IF (RNCYC(I) .LE. 2.0) GO TO 401 
F NCYC( I ) = RNCYC ( I ) - l .0 
GC TO 402 

401 PNCYC(I) = RNCYC(I) - 2.0 

402 KIND = C 
GC TO 415 

410 X 3 = SIGMAX( I ) 

X 4 = SIGMINI I ) 

I ND3 = NSTEP(I) 

I NC4 = I ND3 
KMIN = 1 
KMAX = 0 
K31 = 0 
KIND = 1 

RNCYCU ) = RNCYC ( I ) - 1.0 
KE = 0 

415 CYCNO = A I NT (RNCYC (I >+0.5 ) 

CALL DECIDE(X1,X2,X3,X4,KEY,I , CYCNO , KCY GEN ) 

GC TO 1C00 
20CC LMAX = L 

IF (LRMAX ,LT. 4) GO TO 5000 
IF (NCYNC .EQ. 0) GO TO 5000 
C 
C 

C RANGE PAIR COUNT OF RESIDUE SPECTRUM'S 

C 

C 

NR E S=NR E S+ 1 

CALL CECRES(LRMAX,NCYNO) 

GC TO 2 COO 

500C IF (LRMAX .LE. 1) GG TO 3000 

C 

C 

C COUNT THE LAST RESIDUE SPECTRUM - RANGE-PAIR COUNTING 

C KILL YIELD N ADDITIONAL CYCLES 

C 
C 

K K = 0 

R E $MAX = R ES ( 1) 

RESM IN=PES(1 ) 

I M A X= 1 
I N I N= 1 

CC 500 1=2, LRMAX 

IF (RES(I) .LT. RESMAX) GO TO 450 
RESMAX = RES ( I ) 

IMAX = I 
GC TO 500 

490 IF (RES ( I ) . GT . RESMIN) GO TO 500 
RESMIN=RES( I ) 

ININ = I 
500 CONTINUE 

CALL CY CR E S( R ESMAX, RESMIN, 1.0, I NCEX( IMAX)) 

KK = KK+1 
510 J= I MAX-2 

IF (J -LE. 0) GO TO 550 

CALL CYCRESt RES( J ) ,RES ( J + l ), 1. 0 , INDEX (J ) ) 

KK = KK+1 
IMAX = J 
GC TO 510 

550 J = I MI N + 2 

IF (J .GT. LRMAX) GC TO 575 

CALL CYCRES(RES(J-1) »RES(J) ,1.0, INDEX (J-l) ) 

KK = KK+1 
I M IN = J 
GC TO 550 
575 KM A X = KK 
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INC2 = IND 1 
IN C3 = I NO 1 
INC4 = IND1 
KBIN = 1 
KB AX = 0 
K'l = 0 

IF (RNCYC(I) .LE. 2.0) GO TO 401 
BNCYCI I ) = RNCYC ( I ) - l .0 
GC TO 402 

401 PNCYC(I) = RNCYC ( I ) - 2.0 

402 KIND = C 
GC TO 415 

410 X3 = SIGBAXC I ) 

X 4 = SIGMINI I ) 

I N03 = NSTEP(I) 

I NC4 = I ND3 
KB IN = 1 
KBAX = 0 
K3 1 = 0 
KINO = 1 

RNCYCII ) = RNCYC ( I ) - 1.0 
KE = 0 

415 CYCNO = A I NT ( RNC YC ( I ) +0 .5 ) 

CALL DECIDE (XI ,X2 t X3tX4,KEY , I , CYCNO , KCY GEN ) 

GC TO 1000 
20CC LB A X = L 

IF (LRBAX .LT. 4) GO TO 5000 
IF ( NCYNO .EQ. 0) GO TO 5000 
C 
C 

C RANGE PAIR COUNT OF RESIDUE SPECTRUMS 

C 

C 

NR E S=NR E S+ 1 

CALL DECREStLRMAX, NCYNO) 

GC TO 2 COO 

500C IF (LRMAX .LE. 1) GG TO 3000 
C 

c 

c COUNT THE LAST RESIDUE SPECTRUM - RANGE-PAIR COUNTING 

C KILL YIELD N ADDITIONAL CYCLES 

C 
C 

K K =0 

RE$BAX=RES( 1) 

RESB IN=PES(1 ) 

I B A X= 1 
I B IN = 1 

CC 500 1=2, LRMAX 

IF (P.ES(I) .LT. RESMAX) GO TO 490 
RESMAX = RES ( I ) 

IBAX = I 
C-C TO 500 

490 IF (RES(I) .GT. RESMIN) GG TO 500 
RESBIN=RES( I) 

IBIN = I 
500 CONTINUE 

CALL CYCRESI RESMAX, RESMIN* 1.0, INCEX( I MAX ) ) 

KK = KK+1 
510 J= IBAX-2 

IF <J .LE. 0) GO TO 550 

CALL CYCRES(RES(J)»RES(J+1) ,1.0, INCEX(J ) ) 

KK = KK+1 
IBAX = J 
GC TO 510 

550 J = I MI N + 2 

IF (J .GT. LRMAX) GC TO 575 

CALL CYCRES(REStJ-l) , RES( J ) ,1.0, INDEX(J^i) ) 

KK = KK+1 
IBIN = J 
GC TO 550 
575 KBAX = KK 
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C 

c 

C SC RT THE ANALYSIS SPECTRUM TO PROCUCE THE RANGE-PAIR 

C CCLNTEC SPECTRUM 

C 
C 

300C KP = 0 

CC 605 JJ = 1 ♦ NPKS 
KC = 0 

OC 600 I = l.LMAX 

IF (NNSTEP(I) .NE. JJ ) GO TO 600 

KF=KP+1 

KC=KC+ 1 

NSTEPt KP) =KP 

SIGMAX(KP) = CYCLE(Itl) 

SIGMIN(KP) = CYCLE! 1.2) 

RNCYC(KP) = RNECYC(I) 

IF (KC .LT. 2) GO TO 600 

IF (SIGMAX(KP) .NE. S IGMAX ( KP- 1 ) ) GO TO 6C0 
IF (SIGMIN(KP) .NE. SI GM IN( KP-1 ) ) GO TO 600 
595 KF=KP— 1 

R NCYC ( K P ) = RNCYC(KP) + 1.0 
600 CONTINUE 
605 CONTINUE 
KFMAX=K P 

IF (IPRINT .GE. 2) GO TO 104 
WRITE(6,2010) 

2010 FORMAT (lHl»48X»33hR AN GE PAIR CYCLE COUNTED SPECTRUM//) 
WRITE(6 ,22) 

WRITE (6 ,25) ( NST EP ( I) ,S IGMAX ( I ) , S IGMIN( I ) ,RNCYC( I ) , 1 = 1 
1 , KFNAX) 

102 FCFMAT ( 5X ,3F10.2) 

104 CONTINUE 
RETURN 
END 
C 
C 
C 
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SIBROUT INE DECIDE(X1,X2 ,X3 t X4 , KEY , I , CYCNC , K CYGEN ) 
CCM*ON/MDECl/SIGMAX( 4205 ) , S IGM I N ( 4205 ) 
CCNMCN/MCEC2/NSTEP( 4205 ) , LR , KMAX , KM IN, K 2 l 
CCMMCN/MDECR/RES ( 4205),INDEX( 4 205 ) , IND 1 , I NC2 , IND3 , IN 
1C4,K INC 

CCMNCN/MCYG/CYCLE(2Q0,2),RNECYC( 4205 ) , NNST EP ( 4205) 
CCMMCN/MCGDE/L ,LIND 
C 
C 

C THIS SUBROUTINE CECIDES WHETHER OR NOT ThE VALUES XI, 

C X2, X3 , ANC X4 FROM THE ADJUSTED LOAD SPECTRUf' SATISFY 

C ThE RANGE-PAIR COUNTING CONDITIONS 

C 
C 

KFIRST=0 

IF (K31.NE.0) GO TO 11 

10 IF (X3.LE.X2) GO TO 200 

11 IF (X2.GT.X1) GO TO 210 

IF (X2.LT .X4.GR.X3.GT .XI) GO TO 5C0 

150 IF (X2.GT.X3) GO TO 151 

CALL CYCGEN(X3,X2,l.0,NSTEP(I) ) 

C-C TO 152 

151 CALL CYCGEN(X2,X3,1.0»NSTEP(I) ) 

152 X1=X1 
X2 = X4 

IF ( INC3.NE. IND2) LINC=1 
I NC2= INC4 
KCYGEN = 1 

IF (KEY.NE.O) GO TO 110 
RETURN 

210 IF (X2.GT.X4.0R.X3.LT. XI) GO TO 500 
GC TO 150 
2CC X1=X1 
X2 = X4 
IN C2= IN D4 
KCYGEN = 2 

IF (KEY.EQ.O) RETURN 
CYCNC=CYCN 0-1.0 
GC TO 110 
C 
C 

C ACC XI TO THE RESIDUE SPECTRUM 

C 

C 

500 LR=LR+1 

R E S ( LR ) =X 1 
INCEX(LR) = INDl 
X 1 =X2 
X2 = X3 
X 3 =X4 
INCl=IND2 
I NC2=INC3 
I N C3= I N C4 
KCYGEN = 3 

IF (KEY.NE.O) GO TO 110 
RETURN 

110 GC TO ( 1150,1200, 1500) , KCYGEN 

1150 IF (CYCNO.GT.l.O) GO TO 1151 
IF (CYCNG.LE.O.O) RETURN 

1153 CYCNC = CYCNO - 1.0 
GC TO 1152 

1151 IF (LIND. EG. 1) GO TO 1153 

IF (INC3.NE.IND4) GO TO 1153 
RNECYC(L) = RNECYC(L) + CYCNO - 2.0 
CYCNC = 1.0 

1152 IF (KMAX.NE.l) GO TO ill 
X 2 = SI GM IN( I ) 

I ND3 = NSTEP(I) 

IF (CYCNO. GT. 0.0) GO TO 112 
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KRIN = 1 
KRAX = 0 
KCYGEN = 3 
FETURN 

1200 IF (CYCNO. LE. 0. 0 > RETURN 
CYCNO = CYCNO - 1.0 
X2 = S I GM AX ( I ) 

XA = SI GM IN ( I ) 

KFIRST = 1 
GC TO 113 

111 X2 = SIGMAXUJ 
XA = S I GM IN ( I ) 

IF (KFIRST. NE.O) GO TO 113 
CYCNO = CYCNO - 1.0 
KFIRST = l 

113 I NC3 = NSTEP(I) 

INCA = IND 3 
KRIN = 1 
KRAX = 0 
GC TO 10 

1500 IF (KMAX.NE.O) GO TO 1510 
IF (CYCNO. LE. 0. 0) RETURN 
CYCNO = CYCNO - 1.0 

112 XA = SI GM AX ( I ) 

INCA = NSTEP(I) 

KRAX = 1 

KRIN = C 
GC TO 11 

1510 X A = SIGMIN(I) 

INCA = NSTEP ( I ) 

KRAX * 0 
KRIN = 1 
GC TC 10 
ENC 
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SLERCUTINE CECR5S (LRMAX* NCYNO) 

CCMMON/MCGDE/L,LIND 

CCMMCN/MDECR/RES ( 4205), INDEX! 4205 ) , IND l , I ND2, INC2 , IN 
1C4 ,K IND 



THIS SUBROUTINE DECIDES WHETHER CR NOT THE ELEMENTS CF 
THE RESIDUE SPECTRUM SATISFY THE RANGE P A I R-CCUNT I NG 
CONDITIONS 



K = 0 
NCYNO = 0 
XI = RES(i) 

X 2 = RE S ( 2 ) 

X 2 = RE S 1 3 ) 

X 4 = RE S ( 4 ) 

IND1 = INDEX(L) 

I NC2 = INDEX (2 ) 

IND3 = I NDEX ( 3 ) 

INC4 = INCEXI4) 

J = 4 

10 IF (X2.GT.X1) GO TO 100 

IF (X2.LT.X4.0R.X3.GT.X1) GO TO SCO 

150 IF (X2.GT.X3) GO TO 151 
CALL CYCRE$(X3,X2,1.0, IND3) 

GC TO 152 

151 CALL CYCRES(X2,X3,1.0, IND2 ) 

152 NCYNO = NCYNO + 1 

XI = XI 
X2 = X4 
INC2 = I ND4 

IF (J.EC.LRMAX) GO TO 300 
IF ( ( J+ 1 ) . EQ. LRMAX ) GO TO 315 
X3 = RESU + l) 

X 4 = RES(J+2) 

IND 3 = INDEX ( J+ 1 ) 

I NC4 = INDEX (J +2 ) 

J = J + 2 
GC TO 10 

100 IF (X2.GT.X4.0R.X3.LT. XI) GO TO SCO 
GC TO 150 
500 K = K + 1 

RES(K) = XI 
INDEX(K) = I ND1 
J = J + 1 

IF (J.GT. LRMAX) GO TO 330 

XI = X2 
X2 = X3 

V5 = Y4 

X 4 = RES(J) 

INC1 = I ND2 
I NC2 = I ND3 
INC3 = IND4 
I ND4 = INDEX(J) 

GC TO 1C 
300 K = K + 1 

RES(K) = XI 
R E S ( K + 1 ) = X2 
INCEX(K) = IND1 
INDEX( K+l ) = I ND2 
LRMAX = K+l 
RETURN 

315 K = K+l 

PES(K) = XI 
P E S ( K + i ) = X2 
R ES ( K+2 ) = RES ( J + i ) 

INCEX(K) = IND1 
INCEX( K+l) = I ND2 
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index( k+2 ) = Index ( j + 1 ) 

LRNAX = K+2 
RETURN 

330 K = K+l 

RES(K) = X2 
RES (K+l ) = X3 
R ES ( K+2 ) = X4 
INDEX(K) = IND2 
I NDEX ( K+ 1 ) = IND3 
I NCEX ( K +2 ) = IND4 
LPNAX = K+2 
RETURN 
END 
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SLERCUTINE CYCRE S ( Y1 , Y2 ,C YC PF , NS TE PP ) 
CCMNCN/MCYG/CYCLE<200,2),RNECYC{ 4235 ) , NNST c P ( 4235) 
CCMNGN/yCGDE/L,LINO 



THIS SUBROUTINE GENERATES CYCLES FCR THE ANALYSIS 
SPECTRUM FROM DATA SUPPLIED BY SUBROUTINE CECRES 



L = L + l 

CYCLE ( L » 1 ) = Y1 
CYCLE(L t 2) = Y2 
RNECYC(L) = CYCPF 
NNSTEP(L) = NSTEPP 
RETURN 
END 
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RANDOM I Z E SPECTRUM A LQACS FOR STM4X 



I SK I P= 1 00 
I A = 0 
I E=0 
IC = 0 
IC = 0 
I E =0 
I K = C 
I G=C 
I F = 0 
I 1 = 0 
I J = 0 
1 = 1 

ITCTAL=0 
I X = 9 

700 CALL RAN0U< IX,IY,YFL) 

I X = I Y 

IF (IT0TAL.EQ.4201 ) C-0 TO 790 
IF( (YFL-GE-0. ) .ANO. (YFL-LT..1) ) GO TO 75C 
IF( (YFL .GE..1 ) .ANC. (YFL.LT. .2) ) C-C TO 751 
IF((YFL.GE..2).AND.(YFL.LT..3) ) GO TO 752 
I F ( (YFL.GE..3) .AND. (YFL.LT . .4 ) ) GO TO 753 
I F C (YFL. GE.. 4). AND. (YFL.LT.. 5) ) GC TO 754 
I F < (YFL .GE . . 5 ) .AND. (YFL.LT. .6) ) GC TO 755 
IF< (YFL .GE..6) .AND. (YFL.LT..7) ) GC TO 756 
IF((YFL.GE..7).AN0.(YFL.LT..8) » GC TO 757 
IF( (YFL.GE..8) .AND. (YFL.LT. .9) ) GO TO 758 
1F( (YFL .GE..9). AND. (YFL.LT. 1.) ) GG TO 759 

750 IFUA.EQ.2) GO TO 700 
ITCTAL=ITGTAL+1 

I A=I A+ 1 
RTMAX( I ) = 1 . 25 
BTFIN( I )=.ll 
I TYPE ( I ) = 1 
RNN( I )=1. 

1 = 1 + 1 
GC TO 700 

751 IF(IB.EC.4) GO TO 700 
I TCTAL=ITCTAL+1 
16=18+1 

PTM4X( I ) = L . 1 5 
RTM IN ( I ) = .U 
I TYPE < I ) = 1 
P NM I ) = 1 • 

1 = 1 + 1 
GC TO 700 

752 IF < IC. EC. 15) GO TO 700 
ITCTAL=ITCTAL+1 
IC=IC+1 

RTMAX( I )=1 .05 
RTM IN( I )=. 11 
I TY PE ( I ) = 1 
RNMI » = 1. 

1 = 1 + 1 
GC TO 7C0 

753 IF( ID.EQ.44J GO TO 700 
ITCTAL=ITCTAL+1 

I C = I D+l 
RTMAX ( I ) = . 95 
RT* IN ( I ) = . 11 
I TYPE ( I ) = 1 
PNIS ( I )= 1 . 

1 = 1 + 1 
GC TO 700 

754 I F< IE .EC .136) GO TO 700 
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755 

756 

757 

758 

759 

790 

C 

C 

C 

C 

C 



I7CTAl=ITCTAL+l 
I E= IE+1 
RTMAXII ) = . 85 
RTN IN( I > = .11 
I TV PE ( I ) = 1 
RNNU ) = 1. 

1 = 1 + 1 
GC TO 7 CO 

I F ( IK. EG. 250) GO TO 700 

IT0TAL=ITGTAL+1 

I H = I K+l 

RTMAX( I >=.75 

RTHN(I ) = .ll 

ITYPE ( I )=1 

PNN( I )=1. 

1 = 1 + 1 
GC TO 7C0 

I F ( IG.EC.450 ) GO TO 700 

ITCTAL=ITCTAL+1 

I C- = I G + 1 

RTM AX ( I ) = . 65 

RTF INd )=.ll 

I TYPE { I ) = 1 

PNMI ) = 1 . 

1 = 1 + 1 
GC TO 7C0 

I F ( IH.EC.650) GO TO 700 

ITCTAL=ITCTAL+1 

I F= I H+l 

RTMAXI I ) = . 55 

RTy IN( I > = .11 

ITYPEU ) = 1 

RNM I > = 1 . 

1 = 1 + 1 
GC TO 700 

I F { II .EQ.950) GO TC 700 
ITCTAL=ITQTAL+1 
11=11+1 
R TMAX ( I ) = .A5 
RTR IN ( I ) = . 1 1 
I TYPE ( I )=1 
PNN ( I ) = 1 . 

1 = 1 + 1 
GC TO 700 

IF(IJ.EQ.1700) GO TO 700 

ITGTAL=IT0TAL+1 

I o= I J+ 1 

R TM AX d > = .35 

RTMNII > = .11 

ITYPEU ) = 1 

RNM I )= 1 • 

1=1 + 1 
GC TO 7 CO 
CONTINUE 



RANDOMIZE SPECTRUM A LOACS FOR STM IN 



I A = 0 

i e= o 

IC = 0 
IC = 0 
I E = 0 

i y = o 

I G = 0 
I F = 0 
I I = C 
IJ = C 
I K = 0 
I L = 0 
1 = 1 
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773 

760 

761 

762 

763 

764 

765 

766 

767 

768 



I TCTAL = 0 
I X=583 

CALL RANDU( IX»IYtYFL) 

I X = I Y 

IF(ITOTAL.EQ.lOl) C-0 TO 791 
I F( (YFL.GE.O. ) .AND. (YFL.LT . 

I F ( (YFL.GE.. 07) .AND. (YFL.LT, 
IF( (YFL.GE..14) .AND. (YFL.LT 
I F( (YFL.GE.. 21 ) .AND. (YFL.LT, 
IF ( (YFL .GE..28) . AND. ( YFL . LT , 
IF( (YFL.GE. .35). AND. (YFL.LT, 
I F( (YFL .GE.. 42) .AND. (YFL.LT. 
I F ( ( YFL .GE . .49) .AND. (YFL.LT 
IF ( (YFL .GE. .56) .AND. (YFL.LT , 
I F( (YFL.GE. .63) .AND. ( YFL.LT. 
I F( (YFL.GE.. 70) .AND. (YFL.LT 
IF ( (YFL. GE..77) .AND. (YFL.LT. 
IF ( (YFL .GE.. 84) .AND. (YFL.LT. 
I F ( IA.EG.50) GO TO 773 
ITCTAL=ITGTAL+1 
I A=IA+1 
HCLD(I) =0. 

1=1 + 1 
GC TO 773 

I F ( IB. EG. 20) GO TO 773 
ITCTAL=ITQTAL+1 
ie=IB+l 
HC LC ( I ) =-.04 
1=1 + 1 
GC TO 773 

IF(IC.EC.IO) GO TO 773 
ITCTAL=ITCTAL+1 
I C=IC+1 
FCLC ( I ) =- . 08 
1 = 1 + 1 
GC TO 773 

I F ( ID. EG. 6) GO TO 773 
ITCTAL=ITCTAL+1 
I C= ID+ 1 
H C LC ( I ) =- • 12 
1 = 1 + 1 
GC TO 773 

IF(IE.EG.3) GO TO 773 
ITCTAL=ITCTAL+1 
I E= I E+l 
HCLC( I ) =-.16 
1 = 1 + 1 
GC TO 773 

IF(IM.EG.3) GO TO 773 
ITCTAL=ITCTAL+1 
I N = 1 M+l 
HC LD( I ) =-.20 
1=1 + 1 
GC TO 773 

IF ( IG. EC. 3 ) GO TO 773 
ITCTAL=ITOTAL+l 
I G= I G+l 
HCLD( I ) =-.25 
1 = 1 + 1 
GC TO 773 

I F ( IH.EC.2) GO TO 773 
ITCTAL=ITCTAL+1 
I F = I H+ 1 
HC LC 1 I ) =- .29 
1 = 1 + 1 
GC TO 773 

I F ( II .EC.l) GC TO 773 
ITGTAL= ITCTAL+1 
11=11+1 
HCLC(I) — .33 
1 = 1 + 1 
GC TO 773 



06 ) ) 
.13 ) ) 
. 20 ) ) 
.27 ) ) 
.34) ) 
.41 ) ) 
.48) ) 
.55 ) ) 
.62 ) ) 
. 69 ) ) 
.76 ) ) 
.83 ) ) 
1 . 0 ) ) 



GO TO 760 
GO TO 761 
GO TO 762 
GO TO 763 
GC TO 764 
GO TO 765 
GO TO 766 
GO TO 767 
C-0 TO 768 
GO TO 
GC TO 
GO TO 
GO TO 



769 

770 

771 
773 
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769 IFUJ.EC.l) GO TO 773 
ITCTAL=ITGTAL+1 

I J = I J+l 
HC LD< I ) =-.37 
1 = 1 + 1 
GC TO 773 

770 IF( IK.EQ.l) GO TO 773 
ITCTAL=ITCTAL+1 

I K = I K+l 
HCLC(I)=-.41 
1 = 1 + 1 
GC TO 773 

771 I F ( IL.EQ.l) GC TO 773 
ITCTAL=ITGTAL+1 

1 L = I L + 1 
HC LC ( I ) =— .45 
1 = 1 + 1 
GC TO 773 

791 CONTINUE 
C 

C 

C RANCCMIZE SPECTRUM A STMIN LOADS WITH .11LLC1G.I LCAC 

C TO EQUAL NUMBER CF STMAX LOADS 

C 
C 

I A = 0 
I E = 0 
J = 1 
1 = 1 

ITGTAL=C 
I X=4777 

780 CALL RANDU<IX,IY,YFL) 

I X = IY 

IF< IT0TAL.EQ.4201 > GO TO 792 
IFKYFL.GE.O.J . AND. ( Y FL .LT . .2 ) ) GC TO 781 
IF( (YFL .GE..2) .AND. (YFL.LT.l. ) ) GC TO 782 

781 IF( IA.EQ.L01) GO TO 780 
ITCTAL=ITCTAL+l 

I A= I A+ 1 

R TM IN ( J )=HCLD( I) 

1=1 + 1 
J = J+1 
GC TO 780 

782 IF< IB. EC. 4100) GO TO 780 
ITCTAL=ITGTAL+1 

I e=I5+L 
RTMIN( J ) = .ll 
J=J+1 
GC TO 780 

792 CONTINUE 
C 

C 

C ESTABLISH A GROUND CYCLE EVERY HOUR ( 42 EVENTS ) 

C 

C 

CC 795 J= 1 » 4201 » 42 
RTM I N ( J ) =- .08 
795 CONTINUE 
C 
C 



115 



LIST OF REFERENCES 



1. Coffin, L. F. , Jr., "Low Cycle Fatigue: A Review," 

Applied Materials Research , v. 1, p. 129-141, October 1972. 

2. Vought Aeronautics Division Report No. 2-53420/OR-5562 , 

A7 Wing Fatigue Life , by J. W. Darnell, Jr. and J. W. 

Cherry, p. 1.17, 6 August 1970. 

3. Department of Defense Specification MIL-HDBK-5A, Metallic 
Materials and Elements for Aerospace Vehicles Structures , 
February 1966. 

4. Dowling, N. E., "Fatigue Failure Predictions for Complicated 
Stress -Strain Histories," Journal of Materials, v. 7, 

p. 71-87, March 1972. 

5. Manson, S. S. and Hirschberg, M. H., Fatigue Behavior in 
Strain Cycling in the Low- and Intermediate -Cycle Range , 
paper presented at Sagamore Army Materials Research Con- 
ference, 10th, Raquette Lake, New York, 13 August 1963. 

6. Air Force Flight Dynamics Laboratory Technical Report 
AFFDL-TR-74-23, A User's Manual for the Sequence Account - 
able Fatigue Analysis Computer Program ^ by J . M"! Potter 
and R. A . Noble , May 1974. 

7. Sandor, B. I., Fundamentals of Cyclic Stress and Strain , 
p. 64-84, University of Wisconsin Press, 1972. 

8. Agardograph No. 157, The Accumulation of Fatigue Damage 
in Aircraft Materials and Structures , by J. Schijve, 

January 1972. 

9. Tavernelli, J. F. and Coffin, L. F., Jr., "Experimental 
Support for Generalized Equation Predicting Low Cycle 
Fatigue," Journal of Basic Engineering , v. 84, p. 533- 
541, De ce mb e r 196 2. 

10. van Dijk, G. M. , Statistical Load Data Processing , paper 
presented at International Committee on Aeronautical 
Fatigue (ICAF) Symposium on Advanced Approaches to 
Fatigue Evaluation, 6th, Miami Beach, Florida, 13 May 
1971. 

11. Wetzel, R. M. , "Smooth Specimen Simulation of Fatigue 
Behavior of Notches," Journal of Materials , v. 3, p. 

646-657, September 1968. 



116 



r 



INITIAL DISTRIBUTION LIST 



No. 



1. Defense Documentation Center 
Cameron Station 
Alexandria, Virginia 22314 

2. Library, Code 0142 
Naval Postgraduate School 
Monterey, California 93940 

3. Department Chairman, Code 67 
Department of Aeronautics 
Naval Postgraduate School 
Monterey, California 93940 

4. Professor G. H. Lindsey, Code 67Li 
Department of Aeronautics 

Naval Postgraduate School 
Monterey, California 93940 

5. Lieutenant John Scott Atkinson, Jr., USN 
221 West Street 

Winchester, Virginia 22601 



Copies 

2 

2 

1 

1 

1 



117 



; 63791 



Thes i s 

A877 Atkinson 
c * 1 A study of spectrum 

loading and range-pair 
counting method effects 
on cumulative fatigue 
damage. 



Thes is b 7 ^ 1 

A877 Atkinson 

c.l A study of spectrum 

, oa^ ? ^<7 ra n n«-nn*r 

counting method effects 
on cumulative fatigue 
damage. 



